Electroweak nuclear response at moderate momentum transfer 
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We discuss the convergence of the expansion of the nuclear electroweak current in powers of |k|/M, 
where M is the nucleon mass and k denotes either the momentum transfer or the momentum of the 
struck nucleon. We have computed the electron and neutrino scattering cross sections off uniform 
nuclear matter at equilibrium density using correlated wave functions and the cluster expansion 
formalism. The results of our work suggest that the proposed approach provides accurate estimates 
of the measured electron scattering cross sections. On the other hand, the description of the current 
based on the widely used leading-order approximation does not appear to be adequate, even at 
momentum transfer as low as 300 MeV. 
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I. INTRODUCTION 

Many aspects of the recent neutrino data remain puz- 
zling [l|-[4{ . While unconventional explanations have been 
proposed 043; before advocating exotics, more conven- 
tional effects should be analyzed within realistic models. 

One of the main sources of uncertainty in the inter- 
pretation of the available data is the treatment of nu- 
clear effects, which arc often described using oversimpli- 
fied models of both nuclear dynamics and the reaction 
mechanism. 

Most neutrino experiments take data at beam energy 
E v < 1 GeV, where quasielastic scattering dominates. In 
this channel the elementary neutrino-nucleon cross sec- 
tion is determined by the form factors F^, F^, and Fa, 
where N = p and n for protons and neutrons, respec- 
tively. The vector form factors have been measured to 
high accuracy in electron-proton and electron-deuteron 
scattering experiments, whereas the Q 2 dependence of 
the axial form factor Fa is still controversial, as differ- 
ent collaborations reported significantly different values 
of the axial mass Ma [1, H[ • 

Many studies of neutrino-nucleus scattering are based 
on the impulse approximation (IA), the underlying as- 
sumptions of which are that (i) the momentum transfer 
is large enough for the probe to interact with only one 
nucleon, and (ii) final-state interactions, of both statis- 
tical and dynamical origin, between the struck particle 
and the spectators can be neglected. 

In a previous paper flpj . we investigated the limits of 
applicability of the IA scheme in the case of the linear re- 
sponse of infinite nuclear matter to a scalar probe. The 
results of that analysis suggest that, at momentum trans- 
fer below about twice the Fermi momentum, correspond- 
ing to ^500 MeV at equilibrium density, the IA breaks 
down, mostly due to the oversimplified description of the 
final state. 



In this paper, we develop a formalism suitable to ob- 
tain realistic estimates of the lepton- nucleus cross section 
in the region of moderate momentum transfer. As an il- 
lustrative example, we consider uniform nuclear matter, 
consisting of equal numbers of protons and neutrons, at 
equilibrium density (p = 0.16 fm -3 , corresponding to 
Fermi momentum kp = 1.33 fm -1 ). 



II. ELECTROWEAK CROSS SECTION 



The scattering cross section can be written as 



da cx LauW^ 



(1) 



where L^v is determined by the lepton kinematics and 



w»" = j2( o \ jlJ,f \ n )( n \ jv \°) S ( u 



E -E n ). (2) 



In the above equation, u is the energy transfer, |0) and 
\n) are the initial and final states of the target, and 
J M is the current operator accounting for the lepton- 
nucleus interaction. Within our approach, nuclear matter 
is described using correlated states, including the effects 
of the nonperturbative components of nucleon-nucleon 
(NN) forces fllS- 

We have used Eqs. ([1]) and ^ to obtain the electron 
and neutrino scattering cross sections, including the con- 
tributions arising from one-particle-one-hole correlated 
states. In our calculation, antisymmetrization of the fi- 
nal state and dynamical final-state interactions are both 
taken into account. 

At moderate momentum transfer |q| < 500 MeV, the 
nuclear current operator 



(3) 
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where j M is the nucleon current, is usually treated ex- 
panding in powers of the ratio |k|/M and keeping only 
the leading terms. Here, M is the nucleon mass, whereas 
k denotes either the momentum transfer or the momen- 
tum of the struck nucleon. 
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In this paper, we present the results of a systematic 
analysis of the convergence of this expansion, showing 
that contributions beyond leading order may be impor- 
tant and cannot be neglected. A similar analysis, carried 
out within the context of the mean-field approach us- 
ing the formalism reviewed in Ref. j!3l |. is discussed in 
Ref. 



III. NUCLEAR HAMILTONIAN AND 
NUCLEAR WAVE FUNCTIONS 

The correlated states are defined as [TJ E3 

F\n) 



\n) = 



(4) 



where |n) is a Fermi gas (FG) state and the operator F, 
embodying the correlation structure induced by the NN 
interaction, is written in the form 



N 



F(1,...,N) = S J] hv 



(5) 



j>i=l 



where S is the symmetrization operator, accounting for 
the fact that, in general, [fij,fik] 7*= 0. 

The two-body correlation function features an op- 
eratorial structure reflecting the complexity of the NN 
potential: 



E 

TS 



[hs{nj) + 81, s ftT(rii)Sij] F 



TS- 



(6) 



TS 



In the above equation, = |r,j| with r»j = 
is the operator projecting onto two-nucleon states of total 
spin and isospin S and T, respectively, and 



(7) 



The shapes of the radial functions jTS^Jiy) and /tT (?>/"), 
are determined through functional minimization of the 
expectation value of the nuclear Hamiltonian in the cor- 
related ground state. The assumption that the correla- 
tion operator obtained using this procedure can be used 
to generate the whole basis of correlated states is one 
of the main tenets of correlated-basis function (CBF) 
perturbation theory. It relies on the premise that the 
mean field and the correlation structure induced by the 
nucleon-nucleon interaction are largely decoupled from 
one another. The validity of this assumptions is sup- 
ported by the fact that the results of CBF calculation of 
a variety of nuclear matter properties are in quantitative 
agreement with the results of Monte Carlo calculations, 
as well as with the available empirical data. 

The calculations of matrix elements of both the Hamil- 
tonian and the nuclear current between correlated states 
involve severe difficulties. To overcome this problem, in 
our work, we use the cluster expansion formalism [12j . 
which amounts to writing the matrix element as a sum 



of contributions arising from subsystems involving an in- 
creasing number of nucleons. 

The results discussed in the present paper have been 
obtained at lowest order, i.e. including two-body cluster 
contributions only. Using this approximation appears to 
be justified in the context of a calculation of the nuclear 
matter response at equilibrium density [TH fl(| . 

The Hamiltonian employed to determine the correla- 
tion functions includes the kinetic energy term and the 
Argonne vq NN potential ^7|- The resulting correlation 
functions are then used to construct the effective inter- 
action discussed in Ref. UM needed to obtain the nuclear 



matter single-particle spectrum [16[ . Note that the effec- 
tive interaction of Ref. [l8[ includes the effects of three- 
nucleon forces, which are known to be important. 



IV. CURRENT OPERATOR 

The matrix representation of the weak current is the 
sum of the vector part, 

r$ = 7 M (F 1 + F 2 ) - {p + jfS F2 ■ ( g ) 



2M 



and the axial one, 

r% = y^Fa + 75 



q»F P 

1 = T 15* A T 75 M ■ 



(9) 



In general, the upper Xa an d lower Pauli spinor 4> a de- 
scribing an on-shell Dirac particle of mass M and mo- 
mentum p are related through 

4>a = P ' ^Xa, (10) 

where p is a compact notation for p/(Af + E p ), with 
E p = ^p 2 + M 2 . This relat ion allows one to reduce the 
current according to 

xl, </>l)(r$ + rfi (fy^xK^ + A^, (11) 



with 
V° 



Xp Xpi 



= (F 1 + F 2 ) [l + p' • p + ia ■ (p' x p)] 

F 2 ^±A[l-p'.p-ia-(p'xp)], 



2M 



XpXpi 



= (F 1 + F 2 ) [(p' + iax p') k + {p-iax p) k ] 



- F 



2 ( P + P') fc 
2M 



[1 - p' • p - ia ■ (p' X p)] , 
^=^..(p' + p)-Fp^^..(p'-p) ! 
' " F A [a k +p' k (a ■p)+p k (a- p') - (p' • p) o k ] 



XpXpi 



iF A (p' x p) k - F P 



(P' " P) k 



M 



a - Ip p) 



(12) 
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where = \J (E k + M)/2E k and k = p, p' . 

The standard approach is to expand the above cur- 
rents in powers of momentum over mass, keeping only 
the leading order terms: 

V° = F\ 



2.02 GeV, 15 deg, |q| ~ 537 MeV 



_ , , 7J ^ '[ gx (P'-P)l ^ F i IP + P 



V k = (F 1 +F 2 



2M 



2M 



A o = Fa v(p' + p) Fp? ' 2 - P 2 o ■ (p' - p) ( 13 ) 



2M 



2M 2 



2M 



2M 



Within the considered kinematical region, the terms con- 
taining Fp and Fa in the above equations are of the same 
order, due to the hypothesis of the partially conserved ax- 
ial current relating the pseudoscalar form factor to the 
axial one through Fp = 2A1 2 Fa/ (m 2 + Q 2 ), where 
is the pion mass. 

Note that, in our treatment, off-shell effects arc not 
taken into account in the calculation of the current ma- 
trix elements in Eq. ([2]) , but the single-particle spectrum, 
needed to obtain the energy of the one-particle-one-hole 
states, is modified to take into account interactions ef- 
fects. As a consequence, we can consistently use free- 
space nucleon form factors and the kinematics involves 
no ambiguities. A similar, although oversimplified, pro- 
cedure is used in calculations based on the Fermi gas 
model, in which binding is described through an average 
nucleon separation energy, entering the argument of the 
energy conserving 5 function only. 

In the following, we compare the results obtained at 
leading (LO) and next-to-leading (NLO) order to the ex- 
act result for momentum transfer up to —600 MeV. 



RESULTS 



Figure 1(a) shows a comparison between our results 
and the inclusive electron scattering cross section of nu- 
clear matter at |q| = 537 MeV, obtained in Ref. [19j by 
extrapolating the available data for finite nuclei. The 
curve labeled LO corresponds to nonrelativistic kinemat- 
ics and nonrelativistic single-particle spectrum, while the 
one labeled NLO has been obtained using relativistic 
kinematics and kinetic energies, and including next-to- 
leading-order terms in the expansion of the current. The 
exact calculation includes all terms. Note that most ex- 
isting studies of the nuclear matter response carried out 
within nuclear many-bod y th eory using CBFs are based 
on the LO approximation |20l . [2l| . The results of Fig. 1(a) 
show that, even at momentum transfer as low as —500 
MeV, the contribution of higher-order terms is not neg- 
ligible, and must be taken into account. It clearly ap- 
pears that the NLO and exact calculations accurately 
reproduce the data in the region of the quasielastic peak. 
For comparison, we also show the corresponding results 
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FIG. 1. (Color online) Panel (a): Differential cross section 
da/dudQ for electron scattering off nuclear matter at beam 
energy 2.02 GeV and scattering angle 15 degrees. The exper- 
imental data are taken form Ref. [19j |. Panel (b): Comparison 
of the nuclear matter cross sections with the 2e;Fe(e, e') data 
form Ref. [221 ] . collected at the same kinematics. Calculations 
include the correction for neutron excess in the iron nucleus. 



obtained using the FG model with pure kinetic energy 
spectrum. Note that, while the wrong position of the 
maximum of the FG cross section might be fixed intro- 
ducing an average nucleon separation energy, the failure 
to reproduce height and width of the measured cross sec- 
tion cannot be cured without an ad hoc modification of 
the Fermi momentum. 

To illustrate the importance of finite size effects, in 
Fig. 1(b) we compare our nuclear matter results, in the 



same kinematical setup as in Fig. |l(a)| to the cross sec- 
tion measured using an iron target~[22j . The theoretical 
results have been corrected to take into account neutron 
excess in iron, the influence of which turns out to be neg- 
ligibly small. It clearly appears that finite size effects arc 
large and that a meaningful comparison can only be done 
using extrapolated data. Unfortunately, however, the ex- 
trapolation procedure developed in Ref. [l9[ is based on 
the IA scheme, and cannot be employed for momentum 
transfer below -500 MeV 

In the following, we only discuss the results of nuclear 
matter calculations, for both electron and neutrino scat- 
tering, to illustrate the convergence of the expansion of 
the current operator in different kinematical conditions. 

Figure [5] shows the electron scattering cross sections 
at fixed scattering angle, 30 degrees, and beam energies 
ranging between 600 and 1000 MeV, corresponding to 
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FIG. 2. (Color online) Differential cross section da/dudQ. for electron scattering off nuclear matter calculated exactly (solid 
line), in leading order (dashed line), and in next-to-leading order (dotted line). The panels are labeled with beam energy and 
scattering angle. The values of the momentum transfer refer to the quasielastic peak. 




FIG. 3. (Color online) Comparison of the differential cross section da/dQ 2 for muon neutrino scattering off nuclear matter 
calculated exactly (solid line) and in leading order (dashed line). 



momentum transfer 300 < |q| < 500 MeV. Compari- 
son between the dashed line, corresponding to the LO 
calculation, and the solid line, representing the exact re- 
sult, shows that even at momentum transfer as low as 
300 MeV, the effect of using relativistic kinematics is siz- 
able. Contrary to the nonrclativistic case, when relativis- 
tic kinematics is used, the width of the nuclear matter 
response becomes independent of |q| at large momentum 
transfer, being determined by the Fermi momentum only. 
A suppression of the cross section is also clearly visible 
over the whole energy loss range. 

To make a connection with an observable that is of- 
ten analyzed in neutrino experiments, in Fig. [3] we show 
the Q 2 distribution at fixed beam energies 200 < E v < 
600 MeV. The LO approximation appears to be quite 
accurate at the lowest neutrino energy, where the cross 
section picks up contributions at momentum transfer 
170 < |q| < 300 MeV. As E v increases, inclusion of 
higher-order terms leads to the appearance of sizable cor- 
rections, affecting both the magnitude and the shape of 



the Q 2 distributions. 

VI. CONCLUSIONS 

We have developed a theoretical approach allowing for 
a consistent treatment of electron- and neutrino-nucleus 
interactions in the region of moderate momentum trans- 
fer, in which the IA is expected to break down. Our re- 
sults suggest that the widely used LO approximation for 
the nonrelativistic reduction of the nuclear current con- 
spicuously fails at momenta > 300 MeV. On the other 
hand, inclusion of higher-order terms leads to a remark- 
ably good agreement between theory and extrapolated 
electron scattering data at |q| ~ 500 MeV. 

The main feature of our work is the inclusion of 
nucleon-nucleon correlations, the effects of which are 
long known to be important in electron-nucleus scatter- 
ing [2J]. The use of nonrelativistic correlation functions 
in the final state may be questionable at large momen- 
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turn transfer. However, as |q| increases the role of corre- 
lations between the struck nucleon, carrying a large mo- 
mentum, and the spectators becomes less and less impor- 
tant, and the impulse approximation limits is smoothly 
recovered |10| . 

A systematic comparison with the available data will 
require the extension of our formalism to the case of hnitc 
nuclei. In addition, the contribution of processes involv- 
ing meson-exchange currents, which is known to be sig- 
nificant in the transverse electromagnetic response [25j, 
needs to be included. 

The approach based on correlated basis states, dis- 
cussed in this paper, can be readily generalized to de- 
scribe finite systems. It has been already applied to 
obtain a variety of properties of medium-heavy nuclei 
( 12 C, 16 O, 40 Ca, and 48 Ca), such as ground-state ener- 
gies, charge-density and momentum distributions, natu- 
ral orbits, occupation numbers, quasihole wave functions, 
and spectroscopic factors 12611 . The extension of the tech- 
niques developed in Ref. [26| to the calculation of transi- 
tion matrix elements of the electroweak current does not 
involve any additional conceptual problems. 

The role of meson-exchange currents in neutrino- 



nucleus interactions has long been recognized [27| . Their 
effects must be included in a consistent fashion, as 
the longitudinal components are related to the nu- 
clear Hamiltonian, determining the correlation func- 
tions, through the continuity equation. The formalism 
to obtain the model-independent part of the currents 
within nonrelativistic many-body theory is discussed in 
Ref. [28jj, while the relativistic expressions, as well as 
their reduction, are analyzed in Ref. 13j. It has been 
recently suggested [13, l29U3]| that the effects of meson- 
exchange currents may explain at least part of the dis- 
crepancies between the results of theoretical calculations 
and the double-differential charged-current quasielastic 
cross section measured by the MiniBooNE Collabora- 
tion 0. 
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